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We present the newest statistical and numerical analysis of the matter and cosmic microwave back- 
ground power spectrum with effects of the primordial magnetic field (PMF) included. New limits to 
(— i , the PMF strength and power spectral index are obtained based upon the accumulated data for both 

■ the matter and CMB power spectra on small angular scales. We find that a maximum develops in 

the probability distribution for a magnitude of the PMF of \B\ | = 0.85 ± 1.25(±lcr) nG on a comov- 
, ing scale of at 1 Mpc, corresponding to upper limits of < 2.10nG(68%CL) and < 2.98nG(95%CL). 

While for the power spectral index we find ub = — 2.37~^Q'^(±la) , corresponding to upper limits of 
< — 1.19(68%CL) and < — 0.25(95%CL). This result provides new constraints on models for mag- 
netic field generation and the physics of the early universe. We conclude that future observational 
programs for the CMB and matter power spectrum will likely provide not only upper limits but also 
lower limits to the PMF parameters. 



o 

u 



^ , PACS numbers: 98.62.En,98.70.Vc 

6 

52 ' A. Introduction 

CZ . 

Magnetic fields are everywhere in nature, and they affect many physical processes over a wide range of scales in 
both space and time. Indeed, stochastic magnetic fields have been observed to exist even on the scale of galaxy 
■ clusters [ll-Q and their origin is a recent focus in cosmology fil- flll. Although cluster magnetic fields may be explained 
by a dynamo mechanism, for example, occurring in AGN iets [l6j or via the Weibel instability [l7|], its existence on 
f^*) ■ large scales allows for the possibility that a primordial magnetic field (PMF) could also be present. 

Such a PMF could have influenced a variety of phenomena in the early universe such as the cosmic microwave 
. background [l8l - l3ll |. non-Gaussianity and low multipole anomalies of the CMB[3(], HH, |33| . the formation of large- 
scale structure (LSS) [34rl40| |. and the gravitational wave backgroun d [ill . Effects of a PMF on the BB polarization 
mode of the CMB have also been studied by several authors [ill H3, l23l. l2rl |4^ - I45j and these effects have proven to 
be as important as the other physical processes which strongly affect the BB polarization through Lorentz symmetry 
> I violation HIBD]- 

The purpose of this paper is to present the newest results of constraints on cosmological parameters by taking 
account of the effects of the PMF inferred from the Markov Chain Monte Carlo (MCMC) analysis of the observed 
CMB anisotropies and matter power spectrum. We show that presently existing and future precise observations of 
the CMB anisotropies have the potential to place stringent constraints on the PMF parameters. Indeed, existing data 
may already provide some evidence that a PMF existed in the early universe. If confirmed by future observations, 
this PMF could be a new probe of physics before the CMB epoch and would impact our understanding of the physics 
of the early universe. 

B. Model 

Recently, several different groups (23U27I I29M31I . HH, H(| have been developing the complex analysis of PMF effects 
on the CMB and matter power spectrum. The work presented here, however, is the first complete statistical analysis 
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of the constraints on PMF parameters. Although, occasional discrepancies among different groups have occurred their 
sources have been identified and the groups are now concordant where comparisons can be made. 

To analyze the role of the PMF, one can assume that the photons and baryons behave as a single fluid since 
these particles interact rapidly before the last scattering of photons. Briefly stated, the PMF directly affects ionized 
baryons through the Lorentz force. Since the baryons and photons are tightly coupled before the epoch of last 
scattering, photons are indirectly influenced by the PMF. The cold dark matter (CDM) is also indirectly affected 
by the PMF through gravitational interaction. We assume that the PMF was generated some time during the 
radiation-dominated epoch. In the case of a flat Friedmann-Robertson- Walker (FRW) background cosmology with 
linear perturbations, the time evolution of the energy density of the PMF can be treated as a first order perturbation 
with a stiff source. In this case, all of the back reactions from the fluid can be discarded. For a PMF that is 
statistically homogeneous, isotropic and random, the fluctuation power spectrum can be parametrized [2(| as a power- 
law S(k) = (B(k)B* (k)) oc fc" B where ub is the power-law spectral index of the PMF that can be either negative or 
positive depending upon the physical processes of magnetic field creation. A two-point correlation function for the 
PMF can then be defined[2(| for k < kc where kc is a cutoff wave number in the magnetic power spectrum, i.e. 
(B l {k)B^*(k')) = ((27r)" B+8 /2fc^ B+3 )[B2/r (™B_t3)]fc«Bpy (k)S(k - k'), where P"(Jfe) = 5 ij - ^ . Here, B\ = |B A | 
is the magnetic strength in comoving coordinates derived by smoot hing over a Gaussian sphere of radius A — I Mpc. 
The quantity k\ — 2n/X. Radiative viscosity damps the PMF [49L l50j . We use this effect to evolve the cutoff wave 
number, kc- We then compute the PMF power spectrum using the numerical methods developed in Refs. (26l |38|. 

In this article we place new limits on the magnetic field strength together with other cosmological parameters using 
a Markov chain Monte Carlo (MCMC) method in the numerical analysis jsj. We use adiabatic initial conditions for 
the evolution of the CMB anisotropy and matter power spectrum in the presence of a PMF. Details of the initial 
conditions are summarized in Ref. [23l. |31|. For the scalar mode, we use adiabatic initial conditions for the matter 
contribution as in Ref. [3l| . Although in previous numerical estimations the curvature perturbations of the scalar 
mode were too small to stabilize the numerical calculations for larger scales at early times, our present initial conditions 
stabilize the numerical scheme for all scales and times. Hence, we can now obtain reliable numerical results [5 2]. 

We consider a flat ACDM universe characterized by 8 parameters, i.e. {flbh 2 , f2 c /i 2 , tc, n s , log(10 10 A s ), A t /A s , 
(|i?A|/3.0(nG)) 4 , iib}, where Q,bh? and ft c h 2 are the baryon and CDM densities in units of the critical density, h 
denotes the Hubble parameter in units of 100 km s -1 Mpc -1 , tc is the optical depth for Compton scattering, n s is 
the spectral index of the primordial scalar fluctuations, A s is the scalar amplitude of primordial scalar fluctuations 
and A t is the scalar amplitude of the primordial tensor fluctuations. We define the tensor index of the primordial 
tensor fluctuations as n* = —(A s /A t )/8. For all cosmological parameters we use the same priors as those adopted in 
the WMAP analysis H- 

C. Results 

In our likelihood analysis, the MCMC algorithm was performed until all of the cosmological parameters were well 
converged to the values listed in Table 1. With the inclusion of a PMF, the minimum total x 2 changes from 2803.4 
to 2800.2 corresponding to a change in the x 2 P er degree of freedom from 1.033 to 1.031. Hence, the goodness of 
fit is slightly improved by considering a PMF even after allowing for new degrees of freedom. Figure 1 shows the 
68% and 95% C.L. probability contours in the planes of various cosmological parameters versus the amplitude £>a|, 
or power law index ub, along with the probability distributions. The bottom of the figure shows the probability 
distributions for \B\\ and tib- Of particular note for this article is the presence of maxima for \B\ \ = 0.85 ± 1.25 nG 
and riB — ~2.37^q'y3. Although these values are consistent with zero magnetic field and thus only imply upper limits, 
they suggest the possibility that with forthcoming data (particularly for large CMB multipoles) a finite magnetic field 
may soon be detectable. 

These figures exhibit no degeneracy between the PMF parameters and the standard cosmological parameters. Table 
1 confirms that the standard cosmological parameters are not significantly different from those deduced directly from 
the WMAP 5yr data analysis without a PMF. The reason for this is simple. The standard cosmological parameters 
are mainly constrained by the observed CMB power spectrum for low multipoles I < 1000 (up to the 2nd acoustic 
peak). On the other hand, the PMF dominates for I > 1000. Hence, the PMF effect on the power spectrum is nearly 
independent of the standard cosmological parameters. 

The tensor to scalar ratio A t /A s deduced in our analysis is smaller than the upper limit A t /A s < 0.43 (95% 
CL) deduced from the WMAP 5yr data analysis without a PMF. The reason for this is that we define A t as the 
tensor amplitude of the primary CMB spectrum (without a PMF). This tensor term only arises from the primordial 
gravitational-wave background produced during inflation. We combine the tensor amplitude from the PMF and A t 
when we compare our tensor amplitude with the result deduced by others. The value of A t by itself is comparable 
to the tensor contribution from the PMF. Since the value of A t /A s from the WMAP 5yr data is less than 0.43 (95% 
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CL), our result, is consistent with the previous constraints when the additional PMF contribution is included. 

The degeneracy of PMF parameters [25[ is broken by the different effects of the PMF on both the matter power 
spectrum and the CMB power spectrum. The vector mode can dominate for higher £ in the CMB temperature 
anisotropies [26j . while the matter power spectrum becomes sensitive (38j to the power law spectral index ub when a 
PMF is present. Additionally, the PMF fluctuations are smaller than the CMB fluctuations for the scalar mode on 
larger angular scales. Therefore, the tensor to scalar ratio is not affected by the PMF. 

Figure 2 shows our deduced probability distributions and the la and 2a (68% and 95% C.L.) probability contours 
for the resultant cosmological parameters, as, Hq, z ro i on , and Age. Here, the "alternative normalization parameter" 
<78 is the root-mean-square of the matter density fluctuation within a comoving sphere of radius 8/i -1 Mpc. Hq = 
lOOh km s _1 Mpc -1 is the Hubble parameter, z ro i on is the red shift at which re- ionization occurs, and Age is the 
presently observed age of the universe in Gyr. It is important to keep in mind that these parameters are not input 
parameters, but are output results. The concordance of the results on Figures 1 and 2 provides convincing evidence 
that both upper and lower limits to the parameters of the PMF can be deduced statistically. 

Table 1 summarizes the upper limits to the PMF parameters along with input and output as, Hq, z rc ion, and 
Age. In particular we find that |B A |< 2.10(68%CL) nG and < 2.98(95%CL)nG and n B < -1.19(68%CL) and 
< — 0.25(95%CL) at a present scale of 1 Mpc. Although previous work [25l I2H l38l. |40| could obtain a less stringent 
upper limit to \B\\ they could not constrain tib at all. Moreover, our deduced probability distributions suggest that 
a finite PMF provides the best fit. 

On angular scales smaller than that probed by the CMB the observed number density of galaxies is a better 
measure of the power spectrum. ThereforejSince the PMF mainly influences the small angular scales, using both the 
CMB (WMAP 5yr||, ACBAR 55], CBI 56], Boomerang [H) and the LSS (2dFDR,[58j]) observational data, we can 
constrain the PMF better than in previous works which relied on the CMB data only. In particular, the upper limit 
on «b can be constrained for the first time and the lower limit to ?t-b is approaching the la confidence level. 

The right-bottom panel of FIG. 1 shows that the maximum likelihood is for jib = — 2.2^'g. It is important [25L Elj 
to constrain ub as this parameter provides insight into models for the formation of the PMF. If the PMF were formed 
during inflation one would expect a value of ijb = —3. The value deduced in Fig. I and Table 1, are consistent with an 
inflation generated PMF at the la confidence level. The optimum values of hb, however, if confirmed would suggest 
that the PMF could have been generated at a later phase transition (such as the electroweak transition) or even 
generated after the epoch of BBN [H, El- 

This constraint on ub derives from the gravity waves produced along 
with the PMF. Big-bang nucleosynthesis (BBN) depends upon a balance between the nuclear reaction rates and the 
expansion rate of the universe. Since gravity waves contribute to the total energy density they affect the expansion 
rate. Hence, they are constrained by a comparison between the BBN predictions and the observed light element 
abundances [25|, [4JJ. Forming a PMF at very late times (i.e. after BBN) is a challenge. However, models have been 
proposed (e.g. via second-order MHD effects during recombination [l^]) for the generation of magnetic fields at late 
times which are not excluded from the range of PMF parameters of the 2 a confidence limits deduced here. 

A PMF affects not only the temperature fluctuations, but also the polarization of the CMB. Although we fit all 
available polarization data, it turns out that the TT and BB modes (where T is the temperature fluctuation and B is 
the curl-like component of polarization) are most important. Figure 3 shows a comparison of the computed best- fit 
total power spectrum with the observed CMB spectrum. Plots show various spectra for the TT and BB modes. 
We plot both the best fit and allowed regions both including the SZ effect (scattering from re-ionized electrons) at 
the K(22.8GHz) band (upper curves) and without the SZ effect (lower curves). Including the SZ effect only slightly 
diminishes the best fit magnitude of the PMF. Although the CBI point falls about I a above the best fit, the \ 2 is 
dominated by the better ACBAR08 data and this point does not significantly affect the deduced PMF parameters.. 

D. Summary 

To summarize, this work presents the first complete statistical analysis of constraints on parameters characterizing 
the PMF based upon available data on the CMB and matter power spectrum. For the first time we constrain the 
spectral index ub, which provides new insight into models for the formation of the PMF. The matter power spectrum 
deduced here is consistent with observational constraints on the growth of linear structure as is also reflected in 
the deduced as parameter. We conclude that future plans to observe the CMB anisotropies and polarizations for 
higher multipoles £, e.g. via the Planck, QUIET, and PolarBear missions, we will be able to constrain the PMF more 
accurately. This will permit a better understanding of the generation and evolution of the PMF and provide new 
insight into the formation of LSS as well as a possible new probe of the physics of the early universe. 
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FIG. 1: Probability distributions and contours of la and 2a confidence limits for the standard cosmological parameters as a 
function of the PMF field strength \B\\ and power law index hb- Blue contours show 1 a (68%) confidence limits and sky blue 
contours show 2 a (95%) confidence limits. Red curves in the middle and bottom of the figure show the probability distributions 
of each parameter. Note the existence of a maximum in the probability distributions for |_Ba| and ub- The standard cosmological 
parameters do not have a degeneracy with the PMF parameters because they are mainly constrained by the observed CMB 
data for £ < fOOO (up to the 2nd peak), while the PMF is mainly influenced by the power on smaller angular scales and higher 
multipoles, I > 1000. 
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FIG. 2: Probability distributions and contours of la and 2a confidence limits for the parameters erg, Ho, 2 r eion and Age. Blue 
contours show 1 <r(68%) confidence limits and sky blue contours show 2 <r(95%) confidence limits. Red curves show probability 
distributions for each parameter. Note that these are not input priors, but output results. 



Cosmological Parameters 


Parameter 


mean 


best fit 


n b h z 


0.02320 ± 0.00059 


0.02295 


n c h 2 


0.1094 ±0.0046 


0.1093 


TC 


0.087 ±0.017 


0.082 


n s 


0.977 ±0.016 


0.970 


ln(10 10 A) 


3.07 ± 0.036 


3.06 


A t /A a 


< 0.170(68%CL),< 0.271(95%CL) 


0.0088 


\B x \(nG) 


< 2.10(68%CL),< 2.98(95%CL) 


0.85 


ns 


< -1.19(68%CL),< -0.25(95%CL) 


2.37 


as 


n ci9+ u - u ^ 

U.OlZ_ 033 


0.794 


H 


73.3 ± 2.2 


72.8 


•2reion 


10.9 ± 1.4 


10.5 


Age(Gyr) 


13.57 ±0.12 


13.62 



TABLE I: PMF parameters and ACDM model parameters and 68% confidence intervals (A t /A 3 is a 95% CL) from a fit to the 
WMAP[li] ± ACBARdl ± CBl[5(| ±Boomerang[57|] ± 2dFDR[^] data. 
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